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I. Introduction 


This is the second year of the JPL - sponsored Mathematics Clinic at CGS 
concerned with the metal-oxide-semiconductor-field-effect-transistor (MOSFET)^ 
in addition to the 1984 summer project. The primary focus of last year's 
Clinic was to obtain values of various device parameter constants used in 
one-dimensional (1-D) current models from data supplied by JPL ("parameter 
extraction"), and so the work was numerical in nature. 

The goal of this year's MOSFET clinic was to deliver to JPL a 
mathematical model of the device dynamics from which an accurate and 
computationally efficient drain current expression could be derived for 
subsequent parameter extraction purposes. The initial plan was to obtain a 
simple 2-0 model. However, a careful study of several 1-D models (see [2], 
[3], [4], [7], [11]) revealed many weak points in their derivation. Moreover 
no one of them included an analysis of the source and drain regions. 
Consequently the team decided to look for a more "acceptable" 1-D model, 
namely one whose derivation did not include procedures and approximations 
which we could not justify, and in which the role of the source and drain 
regions would be incorporated . At the same time we wanted to provide a 
comparison of the 1-D models previously mentioned together with some 
explanation of their failure to provide accurate drain current values at 
^ 2 and s -3 (phenomena that were detected by last year's clinic). We 
feel that all these goals have been achieved. 

The 1-D model provided in 3.4 is mathematically sound and 
satisfactory. It allows the source and drain to operate in different regimes; 
it is not based on a questionable derivation of the depth of the depletion 
layer and it includes the contribution of current from the pinched-off part of 
the device, a feature not seen in any of the previous models. 


I 
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The analysis of source and drain regions of 4.2, 4.3, although not 
entirely original is assembled here for the first time. The resistances 
and Rq of the source and drain regions are estimated and used in 4.4 to 
evaluate the potentials and at the source and drain end of the 

channel respectively: 


There are reasonable doubts that R^I and R^I may be related to the 
quasi-Fermi level for electrons at the source and at the drain more than to 
$20 and • Future investigation is needed to decide their proper place 
in the current expression. 

The qual itative and quantitative comparison of 5.1 - 5.5 will provide JPL 
with a useful quick-reference analysis of the most celebrated long-channel 
models, from the Pao-Sah double integral formula to the Brews [2], [3] and 
Wiele [11] closed-form current expression. No extensive introduction to these 
previous models will be made here; reference will be made to the 1984 summer 
report [5] by Morris-Everson and to the comprehensive book by Sze [10]. 

We conclude with the list of symbols and notations (Appendix 1) used in 
the report and a comparative list of them (Appendix 2) as used by the 
referenced authors. 
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II. CHANNEL MODELLING 

2.1 Charge Distribution 

Estimation of the number of electrons in the inversion layer at any 
point along the channel is essential to an accurate evaluation of the current 
at that point. The number may be given by the statistical mechanical distri- 
bution dependent upon potential and thermal energy. One assumes the Fermi- 
level to lie somewhere within the forbidden region, and that there will 
be few enough conduction electrons that the Boltzman approximation to the 
Fermi-distribution may be used. Within the bulk one applies the law of 
mass action to obtain that the product of free carriers will be a constant, 
dependent only upon temperature. That is, 

2 

2.1.1 np = n^ 

Assuming that for p-type material the number of positive ion impurities 
is much greater than negative impurities, N^ » N^, we approximate 

3<l>c 

2.1.2 N^^ n.e , 

and 

2.1.3 NpN^ = n? 

Here N^ and N^ are estimated concentrations. In formulating distributions 
for free holes and electrons, we assume that for a p-type semiconductor, 
the holes (being the majority carrier) feel no effects due to the local quasi- 
Fermi level , so that 

2.1.4 (b — 

Tp Tp 
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and 

2.1.5a 

2.1.5b 


P = n.e 


n = n^.e 


-p((j)-Vg2-<('p) 

P(*-Vbs-*„) 


where all potentials are in volts. 

The inclusion of the body-to-source potential in these distributions 
is to account for the fact that the Fermi -level will be adjusted downwards 
following the application of negative potentials to the body. The inclusion 
of a quasi-Fermi level for electrons allows for variation in electron distri- 
bution along the channel, such that the relative (and complementary) per- 
centages of current due to drift and diffusion may change. The quasi-Fermi 
level is assumed to lie within the interval 


2.1.6 


“^F " '^BS ^ '^n ^ '*’F ■ '^BS '^DS 


The left hand side represents the value of at the source end of the 
channel, while the right hand side represents its value at the drain end. 
Hence 


P(4>-4>r) 

n = n.e at the source 

2.1.7 

n = n.e at the drain 


This allows for variation of the channel potential dependent upon the poten 
tial applied to the drain, The electron concentration is assumed to 

be negligible at any point in the channel when the exponent in (2.1.5b) 
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becomes zero (see 3.4). For every y€ [0,L], where L is the length of 
the channel, we obtain a value x(y) at which <|> = + Vgg. This value 

^ represents the depth of the inversion layer at least in certain regimes 
of operation (see 3.4 for additional details). 

2.2 Poisson's Equation and Surface Potential 

Using the charge distribution of 2.1 we can formulate a partial dif- 
ferential equation satisfied by the electric potential (p. Since we are 
interested only in the stationary state, Maxwel 1 ' s equations are reduced 
to the single Poisson's equation 

2.2.1 div(EgV(|)) = p 

where is the semiconductor permittivity and p is the charge density. 

In our case we have 

2.2.2 p = q(n-p + N' - nJ) 

where and Nq are the ionized acceptors and donors respectively. Using 
(2.1.5a), (2.1.5b) and the assumption that charge neutrality must exist 
in the bulk of the device we obtain 

-p(^-(j)p-Vg5) -3<j)p 

C.C.6 p - -e +e -e ; 

where n^ is the intrinsic carrier concentration. 

If we assume that is position independent (homogeneous material) 

we obtain from (2.2.1) and (2.2.3) 

qn. p((ji-i))^-Vgg) -p(<j)-«(>p-Vgg) p(j)p -3<^p 

(e -e +e -e ) . 


2.2.4 
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If the variables x,y,z are used to denote the depth, length and width 
of the device respectively we have that (p is independent of z and we 
obtain 

^ ‘^'^i . -3(<l>-4>p"Vg2) p(j)p -p<j)p 

— S + — * = (e -e +e -e 1 


2.2.5 


3x^ 3y^ 


) 


We now assume that the horizontal component of the electric field changes 
much more slowly than its vertical component, which implies 


2 . 2.6 


3x2 


at any point of the channel. This assumption (called the gradual channel 
approximation) enables us to write 


2.2.7 


- 5!l(e 

dx2 


”P(<l>“4'p“Vgg) P4*p ”P'l>p 


-e 


This can be integrated simply by noting that ^ 


r r >1 

2,1 

. Idx, 



+e -e ' ) 


2 ^ d_j^ 

dx • 


_d 

dx 


Idx, 


2qn^. ' P(“<l>+<t>p+Vg5) P<l>p "P't’p] 

— -e +e -e 


dx ’ 


2 . 2.8 




r . )2 2qn. 


--X 


Idx. 




e +e +34> 6 ~e 


X=oo 


'I 



p(|.p -pd)pf 


e '* +e +p(t> 

e -e 





_ X 


Our boundary conditions are that (j) ->• Vg^ and d(()/dx ->-0 as x ->• «. 



Simplifying we arrive at 


2 2 9 ^ 

dx 


2qn. 




f 3(0-<D^-Vg5) 3(-<D+(fc+VBc) 


C’VnrJ 3<t>p P'I'p P't'p 

+e r BS ^ p(((,-Vgg)e ^-e )-e 


■3(j)p 3<l>c 


2qn. 






where the dependence on 4ip and Vg^ is implicit. Also, physical considera- 
tions necessitate taking the negative square root as can be seen by the locati 
of d(j)/dx in Figure (2.2.1) . At the surface x = 0 and (() = (Ji^. Therefore 

the electric field at the surface is a function of 4)„. Thus 

s 


dx 


f2qn^l 


x=0 


Pe. 


F(<^g,<^n^ 


It is now possible to establish an implicit 


relation for the surface potential (f)^. Gauss' Law equates the change in 
electric displacement in a region to the charge enclosed in that region: 

J e E • n dA = Z q.. However, we assume there to be no charge present in 


the oxide layer so the field will be constant from gate to surface. The 

electric displacement of the oxide layer, D , can be calculated rather 

ox 

simply and equated with the electric displacement of the body, 


'body 


Re- 


ferring to Figure (2.2.1) we see that (for x-components) 


GS ^s 


"ox -t 


and E 


ox 


body 


2qn. 


P^c 


F(4>s.<}>n) 


Therefore we may equate the two displacements 
D 


- - ) - “01 C F((l) , 0 ) - Di ■ 

ox t ^ GS ^s' ox '^s’^n' body 

ox 


r 


where we have taken C 


ox 


ox t 


and a 


ox 


ox 


2qn,e3 
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Finally 41^ is defined by the implicit expression 


2.2.10 

where 

2 . 2.11 


= Vg 3 - a F((|.3,4.^) 


F(0.s,<t'n) = 




B4>c ”P<)'r $<Pp “P4>„ 


To obtain values for the source end of the channel we use the conditions: 

4>s = <J>2Q and 4)^ = <|ip - Vg^. Similarly, for the drain end of the channel 
we use 4>g = 4>3j^ and 4»^ = 4>p - Vg^ + Vg^. Equation (2.2.11) is the Pao-Sah 
(see 2.4) exact one-dimensional expression for surface potential which can 
be used to evaluate the potentials at the source and drain ends of the channel 
by appropriate choice of 4>^. 

The models of Brews and Van de Wiele presented in the next chapter use 
different forms of (2.2.11) to obtain 413Q and 4»3|^. We shall now describe 
their approach and we postpone to Chapter 5 a comparative numerical study 
of the results obtained. 


Brews 'computation of 4)3^ and 4»3j^ 

In his derivation of surface potential Brews assumes the carrier densi- 
ties (2.1.5a) and (2.1.5b). The surface potential expression itself is 
derived by equating the space charge densities obtained through Gauss' Law 
with the integral of a simplified Poisson's equation (see [3], pg. 19, see 
also pg. 8). That is 
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2.2.12 




(Gauss) , 


2.2.13 


P = -C„a 




+f(*3-Vgs>e 


(Poisson) . 


Equating yields 


2.2.14 


= V 


GS 


P(<l»--ct)„-Vo<;) P<|)p 


1 

2 


P T C 

Since at its largest value, e < 10 ^ << 1, it has been neglected. 

Comparing (2.2.14) with (2.2.11) we see that while the expressions are obviously 

similar there are terms left out. These absent terms are a direct result 

2 2 

of the way Brews takes Poisson's equation as d <j>/dx = q/e^ (n+N^), which 
considers only the effects of acceptor dopant ions and electrons (see [3], 
pg. 5). Accordingly, we find the terms associated with the donor dopant 
ions and the hole density to be missing in (2.2.13). This does not seem 
unreasonable considering their effect; near the surface, (p+Np)<(10 +10 )<<1, whereas 

(n+N^) > (10 + 10 ) . Formula (2.2.14) is clearly just an approximation 
of the exact (2.2.11). It will be seen in Chapter 5 to be an excellent approxi- 
mation "when appropriate". 


Van de Wiele computation of and 

The Van de Wiele model begins a little differently as the carrier densities 
are defined by the following: 


2.2.16 



gP(<>-V(y)) 


and p = N^e 
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P'l’c 

With the relation = n^e we find 


2.2.17 


P(<))-V(y)-(^p) 

n = n^e 


p(-(j)+(j)p+\/gj) 

and p = n^e 


To Van de Wiele, V(y) - Vg^ accounts for the difference between the hole 
and electron qjjasi-Fermi levels ([11], pg. 991). In other words, V(y) - 
Vgs = ‘J'n " <|ip. Substituting this we find the distributions to correspond 
precisely to our previously derived ones ((2.1.5a) and (2.1.5b)). The general 
expression for surface potential is derived in the same manner as Pao-Sah 
and Brews, namely equating charge densities. According to Van de Wiele (see 
[11], (6) and (3)): 

2.2.18 p(y) = -CQj^(Vg 5 - (|>g) (Gauss) , 


2.2.19 


p(y) = -a C 


ox 






) 3 <l»r 

+3(4's-VBs-l/P)e -e 


(Poisson). 


Equating these yields his expression for surface potential: 


2 . 2.20 


"s = “gs ■ “ 


+e 


^<Pc 

+3(<(>s-Vbs-1/P)® 


We see that (2.2.20) differs from (2.2.11) only in the exclusion of the e 

'34>p _c 

term. The fact that e <10 « 1 makes this reasonable and accounts 


for the discrepancy. Comparing (2.2.20) to (2.2.14) we see the Van de Wiele 
model produces a closer mathematical correspondence to the exact solution 
(2.2.11) than the Brews model. However, for the purposes of approximate 


N|M 
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calculations, Van de Wiele suggests distinguishing between strong and weak 
inversion of the surface (see [11], pg. 993). Strong surface inversion is 
associated with the joint condition V^^ » Vj and '1’^ " ^ '^BS ^ '^1’ 

where is the potential at which the transition between inversion and 

depletion regimes occurs (this can be considered a drain saturation potential); 
and Vj, the threshold potential, is the particular gate potential for which 
\l^ becomes equal to the source potential (this is the minimum gate potential 
to produce an inversion regime). For this case, charge density is approximated 
by; 


2 . 2.21 


p = -a C 


OX 




i 


(Poisson) . 


Equating this with the expression from Gauss, (2.2.18), gives an implicit 


relation for surface potential in a strongly inverted device: 


2 . 2.22 


►ss = “gS • “ 






Weak inversion occurs when V«c « Vx or when V-- » V^ and ()>„ 

bb I bS I n 

^BS ^ '^l- 

In this case charge density is approximated by 
2.2.23 p = -a C 


ox 


1 ^*^F 

- ''bs - 


(Poisson) . 


Equating this with (2.2.18) gives an implicit relation for surface potential 
in a weakly inverted device 


2.2.24 


*sw ° “gS ■ “ 


1 ^“^F 

8<*sw - ’'bS - ?)* 


i 


Comparing (2.2.22) and (2.2.24) with (2.2.20) we notice missing terms. 
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Neglecting e -e in (2.2.17) seems quite reasonable since 

-5 

^ they are < 10 « 1. To understand how (2.2.22) and (2.2.24) are derived 

we must consider the charge densities. The total space charge density is 
the sum of the free carrier and depletion charge densities, i.e., = pj^ 

+ Pq. p^j can be considered the effect due to electrons and p^ the effect 
due to holes and doping. We are now approximating the total space charge 
by 

3(*^_“4_”Vgc) 1 i 

2.2.25 = -a C„^(e '►S(V''sS ' P > = "N * '=0 • 

In a weakly inverted device we neglect the effect of free carriers which 
causes the first exponential to drop out, hence (2.2.23) and (2.2.24). In 
a strongly inverted device we assume the free carriers to be much greater 
than the holes and so the second exponential drops out, thus (2.2.21) and 
( 2 . 2 . 22 ). 

In the Brews model, one formula, (2.2.14), suffices for calculation 

of surface potential. In the Van de Wiele model one must first determine 

whether the surface is strongly or weakly inverted to choose the appropriate 

expression. Formulae (2.2.22) and (2.2.24) are clearly approximations of 

(2.2.11) but their accuracy will be seen to be poor relative to Brews. 

To determine the surface potential at the two ends of the channel we 

must consider the change in the quasi-Fermi level, (j>^, based on its position. 

It is assumed that * at the source end of channel and at the source itself 

n 

are equal (see [3], pg. 20). Thus to calculate surface potential at the source 
end of the channel we may replace by (p^ - Vg^ ([3], (4.22)). Simi- 
larly, at the drain end equals at the drain itself and we may replace 
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by (j(p - Vgg + Vpg ([3], (4.24)). These substitutions may be made in 
all three models. To sum up then, the three models generate four similar 
expressions for surface potential at the source and drain ends of the channel. 
At the source end ((J>^ = (Jip - “bs'^ 


Pao-Sah 


2.2.26 


♦so ' “gS ■ “ 


6<*S0'*f) ^(-♦sO**F*''bs) 
e +e 


p(J) -Pljlp 


-b(*f-“bs*''ds> 

-e -e 


Brews 


2.2.27 


■^SO ^ '^GS ‘ “ 




+p(^so"''bs^® 


^<t>t 


Van de Wiele 


*‘SW0 ^ '^GS ■ 


1 

l^^^swo * ''bs - p)® 


2.2.28 


'^SSO ■ '^GS ■ “I® J • 

At the drain end ((|)^ = 4>p - Vg^ 


Pao-Sah 


2.2.29 


^SL " '^GS " “ 




p(j)p -g(})p p(j)p "P(4 >p"V35+Vq3) 

(e -e )-e -e 


i 
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Brews 


2.2.30 


♦SL = ''gS ■ “ 




+8(*sl"''bS*® 


He 


Van de Wile 


“^SWL '^GS ■ “ 




1 


2.2.31 


'*’SSL ' '^GS " 




2 . 3 Derivation of J ^ ^(X3Y) 

The current in a MOSFET is the result of electron flow in the inversion 
layer, which is governed by the drain-to-source potential. This causes a 
drift towards the drain and a statistical diffusion due to heat energy. These 
two effects are summed up as follows: 


2.3.1 


I = 


^^drift -^diff.^^^ 


2.3.2 

2.3.3 


where D 


_ - 3(f) 

■^drift ~ ^^n 3y ’ 

■? - _«n In 

'^diff. “ 9y ’ 

is the diffusion rate determined by the Einstein relation as 


2.3.4 D = # . 

p 

The current density due to drift is proportional to the product of mobility 
of electrons, number of electrons, and electric field along the channel. 
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Because no current leakage is assumed through the oxide and depletion layer 
except for negligible recombination effects, one may assume that the contri- 
bution' of the vertical field component to the current is negligible with 
respect to the contribution of the longitudinal component. 

The current density due to diffusion is proportional to the electron 
concentration gradient 9n/9y where the number n of electrons is deter- 
mined by the potential (|», according to the Maxwel 1-Boltzman approximation 
to the Fermi -Dirac distribution. Here we allow 


2.3.5 


on - on ^ (n 

-qO - -qO (n,e ) 


Hence 


2.3.6 




i = -nn/M - — 

"^diff. *^^^9y 9y ^ 


Equation (2.3.6) provides a much needed relation between the drift and 
diffusion components since both are now expressed in terms of the surface 
potential, <j>, and the electron quasi-Fermi level 4>^. Note that the variation 
in concentration along the channel is accounted for by variations in surface 
potential, «j»g, and quasi-Fermi level, such that the electron distribution 
changes as one approaches the drain. Upon substitution of (2.3.6) and (2.3.2) 
into (2.3.1) we obtain 

2.3.7 I = wjciun dx 

where the integral is over the inversion layer. 
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Knowing the electron concentration, n, and the quasi-Fermi level gradient 
at any cross-section of the channel we can derive the current flow 
^ through the section by integrating (2,3.7). The principle of conservation 
of charge will then allow us to eliminate the y-variable with an integration 
along the channel to obtain 

W ff 

2.3.8 ^ = r JJ ^ ’ 

It is worthwhile to note that since the total current flow at any cross-section 
of the channel is the same, the diffusion portion of the current increases 
as the drift component diminishes and vice versa. The use of quasi-Fermi 
level for the electrons carries this information, by assuming that changes 
with y, while the Fermi level remains fairly constant. 

2.4 The Pao-Sah Formula 

The Pao-Sah formula is an "exact" one-dimensional expression for the 
drain current based on (2.3.8). We call it "exact" in the sense that no 
other simplifications are introduced in addition to the gradual channel 
approximati on. 

It is based upon the idea that at any point along the channel the current, 

I, is equal to the integral of the current density over the inversion layer: 

Xi(y) 

2.4.1 I(y) = W J„(x,y)dx . 

Jq ^ 

Using the expression for current density, (2.3.7), one gets 

(0 di (y) 

2.4.2 I(y) = Wqy n(x,y) dx . 

^x.(y) 
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V 


Since the quasi -Fermi potential for electrons is assumed to be independent 
of X, the y-dependence may be removed by integrating both sides of (2.4,2) 
along the channel. Hence 


2.4.3a 


_ WSM 
■ L 


rL 


d(j)^(y) 

dy 


ro 


x.j(y) 


n(x,y)dx dy 


2.4.3b 



‘*’f''^bs 


fO 

' x^(y) 


n(x,y)dx dc|i^ . 


The limits of the first integral in (2.4.3b) are the value of at the 
source and drain ends of the channel, as postulated in (2.1.4). 

The inner integral of (2.4.3) must now be put in terms of 4)^ and tp. 

This integral represents the charge present in the inversion layer at an 
arbitrary point in the channel, which is known as a function of potential. 

Using the Poisson equation (2.2.7), we may integrate once, using the boundary 
conditions that 4> - and d(j)/dx are equal to 0 at x = ®. This obtains 
the function in equation (2.2.11), or 


2.4.4 


dx 


”F ( 4> , (|)^ , ) 


Using thisfact, we may express the inner integral of (2.4.3) as 


2.4.5 


rO 


x^(y) 


n(x,y)dx = 


4> 

rs 




n((|),4>p,ct>^)dij 

^dx^ 


r s 






^B(rV*F"''BS> 

’ F(»,»^,»|.) ='♦ 


where ('^’p'^Vgg) and 41^ represent the value of 4i at the bottom and top 
of the inversion layer, respectively. 
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2.4.6 


Combining (2.4.5) with (2.4.3), we obtain the double integral formula 
3LpWqyn. j-^s 

I — I f- / ± A A \ 




'^F^^'^BS 




This formulation is general, and valid in all regions of operation. 
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III. SIMPLIFIED LONG CHANNEL MODELS 


3.1 The Pierret-Shields Single Integral Formula 
^ The Pao-Sah double integral formula given in (2.4.6) is; 


(3.1.1) = 




^ — 


’^‘-•'-D <5^(0) "F 


where 


%(0) = - V 


BS 


and 


%(l) = Vos ^ 


BS 


Note that in this section, $ is measured from the substrate. 

The $ integral is essentially over the inversion layer. Since -t is 
small in the depletion region outside of inversion, we can subtract the $ = 0 
value of electron concentration and extend the integration through the whole 
depletion region. Thus (3.1.1) becomes 


(3.1.2) In = 


Wu, 


/2L.L, 


$^(0) '0 F($,$^,$p) n 


with 

(3.1.3) 


0($-<f ) 3($p'<6) 6<tp e«r 1 

F($,$^,$p) = {e + e + 6$(e )-e -e 


The key thing to realise now is that 


iE = fe 


g($.$ ) -g$ 

-e ) 


2F($,$^,$p) 


(3.1.4) 
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Thus (3.1.2) will become an integral with respect to F instead of $ 


(3.1.5) 


VIWu, f^s 

^0 el.Ln , n 

D ‘’$„(0) ■'0 
n ' ' 


At this stage, referring to Figure 3.1.1, we change the order of 
integration to give: 


(3.1.6) 


^ (L) 
■ n 


i„(0) 






4.„(0) 




d<J> 

^n " 


where the curve in Figure 3.1.1 is given parametrically by: 


<l> = $„(y) 


$ = i^(y) . 


For a given "** 5 (y) , specified by the Pao-Sah gate voltage relation 


( 2 . 2 . 10 ) 


22 



Hence we arrive at: 


4 >n(L) 4 * 

r <> r S 


<i>, 


^(0) J 


SL 




SO 


a4>n n 




F(4>,0,4>p)d$ 


Jo 


F(4>g,‘}'n,‘J-p)d-i>s . 

We can evaluate the last term of (3.1.7) by referring to (2.2.10). 
(''gB - »s) = F(*s.»*. »f) . 

''gb = ''gs - ''bs- 

The last integral thus becomes: 
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(3.1.9) — ^ -- ^ d<i>^ 


P*-D^ox rr, /. .X 1 /. 2 .2, 

^2e. '^GB^‘^SL ’ ‘^SO^ ' 2 ^‘^SL " ^SO ^ 


Thus we derived the single integral form of the Pao-Sah current expression: 


W(iC 

(3.1., 0) 




K-*so') 


W(i e .\/F‘ 


F('i^,+^(0),4>p)d^ - F('i>.4>„(L),4>p) d$ 


where <!>_- and <!>(.. are calculated 


(3.1.11) V 




(3.1.12) V5b = «.5^*^ • 
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3.2 Brews Model 

The Brews charge-sheet model is derived with the same reasoning which 
leads to the double integral form mentioned in (2.4). The number of electrons 
in the inversion layer is estimated using the charge-sheet approximation 
and the current density at any point of the channel is found to be 


3.2.1 


I = PQ 


d(p 
^n 

n 3y 


where is the charge of electrons/unit area and <|»^ is the quasi -Fermi 
level for electrons. We now illustrate the method used by Brews for computing 
and 3<j>^/3y. 

To do these calculations we need to know how charge is assumed to be 
distributed within the device; this is virtually the defining characteristic 
of any simplified one-dimensional model. 

The important thing to realize in Brews' derivation, is that he only 
considers two types of charge in the MOSFET; namely free electrons, n, and 
ionized acceptors, N^. 

The omission of ionized donor charge, N^, is acceptable, since it always 
occurs with in the form - Njj. Since Nq « N^, the approximation 
- Np s is seen to be good. 

The omission of charge due to majority carriers (for n-channel devices 
these are holes, p) is more noteworthy. In fact, p is very small at the 
surface, and remains small up to some considerable depth — this is essentially 
the depletion region. Proceeding deeper into the device (x increasing) 
p tends asymptotically to N^, its value in the bulk. Thus there is no 
definite lower boundary to the depletion region. The approximation Brews 
uses is that p = 0 down to some depth x = w, where it attains its asymp- 
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totic value. Thus we can think of the charge being distributed in the following 
manner. At the surface is an infinitely thin "charge sheet" of electrons; 

^ below this there is solely the constant doping charge, N^, down to some depth 
w; beyond this, charge neutrality exists. 

To find the current, we will first deal with Q^. Since direct evaluation 
is very difficult Brews uses the equality 

3.2.2 Q„ = 5s - 5o 

where is the total charge/unit area and Qq is the charge/unit area 
in the depletion layer due to ionized impurities. This is found to be 


3.2.3 


Qn = "PNflW 


with w equal to the depth of the depletion region. At this point we encounter 
one of the most questionable steps of Brews' derivation. To obtain w he 
writes 


3.2.4 


d^d) _ 
dx2 ■ " 


0 < X < w 


which holds in the depletion region before inversion takes place. Therefore, 
using the conditions (p = Vg^, d(j>/dx = 0 at x = w we have 


3.2.5 


qN. 

<J'(x) = Vgg + ^ (x-w) 


Now it is assumed that (3.2.5) remains true after inversion and can be extended 
to find (j>g, the value of iji at x = 0, This gives us 


w = /2 - fVgj) 


i 


3.2.6 
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Brews then states (but does not use) the result 

V 3.2.7 w = n LgO<).g - pVgg - 1)^ 

which he claims (without proof) comes from a "more accurate derivation that 
includes majority carriers more carefully". 

(3.2.3) and (3.2.6) give us 

3.2.8 Qg = -yf qN^Lg(3«.^ - . 


To complete the expression for Q^, we need to find Q^. Continuity 
of the field across the oxide-silicon interface gives us the following: 


3.2.9 


d$ 

'ox dx 


- e 


x=0 


s dx 


x=0 


The left-hand side of (3.2.9) is given by 

e. 


3.2.10 


do 

'ox dx 


x=o- t-S 


Finally, Gauss' Law gives us to be 


3.2.11 

Thus we have 


Q = £ ^ 
^s s dx 


x=0 


* = -':ox(''gS - 


3.2.12 = ^2 qN^Lg(M^ - 


Now we turn to the calculation of d<p^/dy. 

To find this, we relate to using once again the continuity 
of field at the surface (x = 0). The right-hand side of (3.2.9) is found 
by integrating the following initial value problem: 
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3.2.13 


dx2 


= 


+ e 

s 


qn, P(*-Vg3-*„) 


0 ^ X < w 


w < X 


^ = 0 and = V 03 at x = w. 

Evaluation at the surface, x = 0, (ji = (j>^ gives: 


3.2.14 


dx 


= 

xso"^ P>-B L 




M2 p(^3-Vbs) ■ 


With (3.2.4) and (3.2.10), (3.2.9) becomes 



v ^2 e 



2.15 

<^ox'V*s)= SLj 

s<V''bs>* 

[tJ 


2 3(4>c“<<>f,) P(‘i’e"'^Rc) 

e " (e ® - 1)1 




which, after some manipulations gives 
3.2.16 


n h ^n.'i, 


N. 


2 PC-f.-Vn.) 

(e ^ - 1 ) 


with 

3.2.17 


/2 


a = 


^ox'-B 


With 3<|)^/ay obtained from (3.2.15) and given by (3.2.12) Brews 
writes (3.2.1) as 

.2 


3.2.18 


I 


WuC 


ox 


3 






s “BS- 




s ^"BS- 


W 



Rather than integrate this exactly with respect to y, Brews notes that 


the second bracketed term (the fraction) will be important only when the 


first becomes negligible. Since the first represents Q^, it will be small 


only near pinch-off. But, in this case we can write 

2g(Vp--(J) )+a^ 

3.2.19 r s 1 + 


P(Vgs-*s)+a(34>s-3VBs)‘ 


2(3VPVbs) 




Then 


3.2.20 


I = 


WuC 


ox d 


3 dy 


3Vgs-iMo - ^ + 1 + a(34»,-3VBJ" 


's ‘^’BS 


s ^“BS^ 


Integrating (3.2.20) from 0 to L one obtains the source-to-drain 
current in terms of 4>gQ and the potentials at the source and drain 
end of the channel as given by (2.2.27) and (2.2.30) 

WyC, 


3.2.21 


I = 


ox 


L3' 


( i+sVjs ) < s*sl-Mso [ ( e*si_ )^- ( e«so )^] 




+ a[(34>2|^“3Vg2) ^ ^ 


SO ^“BS^ 
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3.3 The Van de Wiele Model [11] 

Van de Wiele uses a quasi-Fermi level formulation for carrier densities 
as in (2.1). 

From (2.2) we see that his general surface potential equation is given 


by 

3.3.1 
with 

3.3.2 


►s = “gS 






F(q3,*„) 


e +e 


POp -P«l»n 


+[P(V'^BS^“^^® -e 


Differentiating (3.3.1) with respect to (j)^ yields 


3.3.3 


do>, 

d(t> 


-P(V'I»F^ 


1 + 


+23K(VQ3-<0g)+l-e 


P(Vbs-^s) 


2 2 

VBs-S*3-e 


where K = 


ox 


Now, p(<|)^ > -p<))p > 11 and p(|»p+P<))j^ ^ P<j>p > 11. Thus e 

. P^''bs"'*’s^ ,.-5 

and e < 10 . 


■^(V'^’f) ^ ^Q-io 


Hence the following is a good approximation to (3.3.3); 


3.3.4 




1 

where M = ^ = — — 
^ox 
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For the current, the same approach as (2.3) is taken, but with an 
effective mobility introduced: 

f'- 4c = 

J n 


3.3.5 


Hence 


3.3.6 


"eff" 


where (|)^(0) = <l>p-Vg 3 and (J>^(L) = 4»p+Vg5-Vg2 as in (2.2) 
Now we use 


3.3.7 


SC ■ % 


with 


3.3.8 


Qsc ~ ^ (field at the surface) = ~ •J'g) 


and 

3.3.9 


Qd = - 






0(V''bs)-i^ 


(3.3.9) is equivalent to (3.2.8) in Brews' calculations. 
Thus (3.3.6) becomes 


3.3.10 




<=ox<V*s> - 


2S->V‘ 


(p<V’'bs>-i> 


d(|). 


By using (3.3.4), we find that: 


31 




32 


Where 


^ ■ '*’SL " '^BS " p 


® ' *^S0 ' '^BS ’ p 


= ''gs - ''bs - 1 * r 


a = 




ox 


The part in the first parentheses is equivalent to (3.2.21) of Brews. 
Numerically, the second term in (3.3.13) is seen to be small compared to the 
first. If the approximation of [11], eqn. (3.2) is used in (3.3.11), it is 
seen that only the first term appears in the final current equation. Hence 
the Brews charge sheet result is an approximation of (3.3.13). 

The derivation here of MOSFET current still depends crucially on the 
result Qr, = as found in Brews (3.2.8). This is the step 

which we feel has not been sufficiently justified. 

The surface potential schemes presented in (2.3) are seen to be inaccu 
rate in some bias regimes. This is not a problem, however, since we can use 
the full one-dimensional equation to find and (see (2.3)). 

Thus we conclude that, with the exception of using the result (3.3.9) 
for Qq, the Van de Wiele derivation of drain current is based on assumptions 
and approximations which are nowhere too harsh. If the result (3.3.9) for 
Qjj were confirmed, then this model would be virtually optimal as a closed 
form of the one-dimensional Pao-Sah model. 
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3.4 The CT inic's Model 


In the following section, we will give the details of a model derived by 
the Clinic which gives the drain current in all regimes of operation. It has 
K not been possible to do a complete test on this model since the approach was 
devised very late in the semester. Some preliminary numerics do, however, 
indicate that this is certainly a valid model. 

The strength of this model is in that, unlike previous models it does 
not use the same approximations over the whole range of operation of the 
device (a formidable task), or even, in fact over the whole device. 

We have divided the device into regions A, B, and C . A full 
description of these regions is now given, together with various 
approximations valid therein. 


REGION A 


This is defined by 

(3.4.1) *s ^ ''bS *F * 

In this region, the surface is strongly inverted. 

Let's look at the relative sizes of charge densities in the region. 


(3.4.2) 


n = n^. e 




(3.4.3) 




= n. e 


(3.4.4) 


p = n. e 


e(*F*V35-4) 


Since we make comparisons with N^, and Nq « N^, we need not consider 

5 

Note that and e > 10 . 
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At the surface, n > n. e , and thus n/N. > 1. Also we have p < n. 

■«»F -10 

e , and thus p/N^ < 10 ! Hence n/N^-p = 


Now consider a point at depth x. beneath the surface where $ falls to 

-S$ -S$p 

the value Vg^ + . At this point p/N^ = e and n/N^ = e ; hence 

-5 

both are less than 10 . (Still we have n/N^^-p = n/N^.) This should now 

suggest the basis of the new approach. 

Seeing how dramatically n/N^-p falls off between x=0 and x^ , we are 

justified in assuming that the contribution of charge due to electrons comes 

almost exclusively from the region 0 s x < x^. . 

The charge in the device below x^. is given by the field at x^ , namely 


. d4> 


2qn^ 

3 dx 

X=X. = 1 

3 J 


In the function F , the significant terms are: 


B(*-VBs-*n) 


6(VVbS-*J 


(g«-BVg5-l)e 


0*r 


At X = x^. , we have that: 

‘nl , ») «(‘F-*n) 

e = 1 : e = e 


and 


S$r B$p 6*r 

(S($-Vg5)-l)e '■ = (e*^-l)e i (B4p-l)e > 10 


Thus at X = x^ 


B*c/2 

(3.4.5) F == e '■ 
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In fact, neglecting the first term amounts to neglecting the contribution 
due to electrons. 

Thus we have 


(3.4.6) 


Qq(A) = - 


B 




5 


the charge beneath x^. due to We will refer to N^-p as depletion 

charge. 

Finally, we use the fact that the potential drops off rapidly as x 
increases from zero. This allows us to neglect the contribution of depletion 
charge from the region 0 s x < x^ . 

Hence Qq is a good approximation for the total depletion charge in the 
device, and using Gauss' Law we find , the charge due to free electrons: 


(3.4.7) Q^(A) = - Qq(A) 

Now we use Qjj. = £5 ^|_(x=0) = "Cox(Vgs”^s) obtain 


(3.4.8) 


Qn(A) - 




ox^’'GS 


We now make some definitions to be used throughout this chapter. 
Definition 3.4.1 

At any point along the channel, the region osx<x^ is defined to be the 
inversion layer. 

Definition 3.4.2 

At any point along the channel, x = is defined to be the depth at 
which = Vgj + . 

Now let us move on to 
Region B 

% " *s " ''bS ^ ^ 


(3.4.9) 
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Once again, we approximate the depletion charge beneath by 

(3.4,6). 

Now, however, the inversion is not as strong. To account for the 
^ possibility of significant depletion charge in the inversion layer, we do the 
following calculation. 

- 3 ($-Vg£) -Btp 

We still have p/N^ = e < e < 10 ^ hence if we can find 

x^. ,then the depletion charge in the inversion layer is just • 


(3.4.10) 


f ^ 

r^n *BS 



dx ■ 

d4> , 

[K‘ 


dF 
. dx 
‘*'s 

2qn. 


.0 j 


4 


*_+V, 


d4> 


(3.4.11) 






2qn. 


A 


d4> 




The significant terms of F are. 


b(^-Vbs-%) e($p+Vg3-<&) 

(3.4.12) e ; e ^ 


once again: 

3*p 

; (3$-3Vg3-l) e 


$ +Vr,c<* ^ ^ we find that: 

n BS s n F B5 


For 
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^ ^ ^b(vVbs-*) 

e$P 6$P 

(0$-SVg3-l)e > (6$p-l)e > lie 


< 1 ; and 


e4*r 


If we retain only the third term, then we find that in the worst 
possible case = *p"Vb 5 ) sul at the worst end of the integral, we have: 


1 1 

F(approx.) ~ F(Correct) ,, q 

F(Correct) 

Thus the approximation of keeping only term #3 is not at all too harsh. 
Hence the extra depletion charge is: 




(3.4.13) 




2qN, 


-U 


(pI'-pVgj-l) * d+ 


V''bs 


(3.4.14) -N^qx^ = -2 { Vp4— pV^j-l - \/p4.^-l) , 


where 2 

Thus in region 
(3.4.15) 






L ^ J 

B, we approximate 
== -Cox(VgS-^s) " 


Qn by: 


Q/B(#s-Vbs)-1 


which is precisely the form used in Brews and Van der Wiele. 

NOTE: We have derived this as a result of approximations which are only 
valid in region B . The approximation for x.. fails totally in region A. 


Finally we look at 
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Region C 
(3.4.16) 


'bs " % 


This is the pinch-off region. The potential has fallen below the effective 
Fermi-level, and thus the electron concentration is very small. 

In this case, it is possible to evaluate directly. 


* 


(3.4.17) 

Q„(C) = - q 

n(x,y)dy = -qn. 

2qn^. 

*5 


4 

0 




e 


4>P+V 


BS 




The lower limit has been set to 


but this is subject to 


debate. 

Referring to (3.4.12) we see that since ♦p * Vjs < ♦ s .Vjs , 

the first and second terms of F are both less than one, whereas the third is 
at least 10^. Thus we neglect all the terms of F except the third. 

Then becomes 




(3.4.18) Q^(c) * 


Pqn.Bj 


-iH, 


2e 








This approximation for F is still good if the lower limit of 

integration is extended to + Vg^. 

2 

Under the substitution t = B($-Vgj)-1, this becomes: 
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(3.4.19) Q^(C) = 


2qriie, 


-i»s 




Pe 




t2 

dt 




Since the equation for surface potential. 


(3.4.20) 


''gs-*s 


S»r -9«. ! 

♦ (6»5-6Vj5-l)e B -e , 


becomes essentially independent of , and thus will be constant in the 
whole pinched-off region, and the numerical integration in (3.4.19) need 
only be done once. 

This concludes the discussion on Regions A, B, and C, in which we found 
simple expressions for Q^. How do we find the current? Most generally, 
under a given set of applied voltages, parts of the device will be in each of 
the 3 regions A, B and C. 


The current is given by 




(3.4.21) 


I = 


L 


Q d# 

n 




Since we have three different expressions for , it will be necessary 
to split the integral into three. 

Consider a typical case where the source is in A , the drain in C , 
and a portion of the device in B . 
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* * 

The boundary between regions A and B is at ♦=♦,<11 = ♦ . This is 

n n $ s 

defined to be the point along the channel where becomes equal to 


Thus to find we must solve the implicit equation (2.2.7) 


♦n * ''bS (3.4.9)) 

* 

Thus to find ♦ we 
n 

★ ★ 

with Namely 

(3.4.22) (V^c-4.% -V 1 = 

'^GS ''bS^ 


2qnie^ 


-ih 


t; '^'*n^F"''8S-<) • 


The boundary between regions B and C is defined to be where 
*s = *n '^BS ■ drain is pinched off, the surface potential in the 

pinched-off region is constant and thus equal to . Hence this boundary 

IS at = ♦SL"'^BS * 

Thus the current is given by: 


(3.4.23) 



^ n 

'bs 

r n ^ 



Q„(B)dk„ . 

0„(C)d»„ 

k 

4 

♦„(0) 

* 

n J 

♦sl-''bs ^ 


Now we do the integrals. 


(3.4.24) 


(3.4.25) 


* 

K 

r n 


Q„(A)d*„ 


* 

f n 


#„( 0 ) 


4>. 


- Cjj^(Vqs- 4>5)] d4>, 


or 


♦„{ 0 ) 


4>. 


Q„(A)d4>n = 2(Vp4>*-l - Vp4>^{0)-1) - C 


ox 


♦„( 0 ) 


4> 


d# 

so 


41 


To evaluate the second term of (3.4.25), we use the approximation for 
^*^n used by Van der Wiele (3.3.4) . The integral thus yields: 

d*s * 

t n 


(3.4.26) 


Q^(A)d<^^ - 2(Vp^ - \/p4>^(0)-l) 


4>^(0) 


P 


ox 

2 


(pV|.5.2)04.*-p4.5o) * \ 




P^<V''gS>^ ■ "'K • 




"SO '"'BS 


a 


a4M.;.pVGs-f , >-P^'so+P«Gs"f 


"-PVP''gs*I “*P*so-^''gs- I 


where M 


2pqN.e 


A'-s 


ox 


rM 


and a « M(^ + ^ 


The integral for region B is precisely that derived in Van der Wiele 
(3.3.13): 


(3.4.27) 


r*SL“^BS 


Qn(B)d4.^ 


4>_ 


4>, 


SL 


d4> 




where 
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(3.4.28) 




Q ^( B ) d.^n 

•^4>* 


OX 


'P^gs-^X-Sl-^so' - f(^s^♦so> 


+ 2 '•/2 1/2 
+ 3 a P(A - B ) + 4 a (A - B ) 


+ In 


1/2 

Vqs 


♦s0-“ Vg3 


where 


1/2 

+ 2 a C fn 


* "SL - ''es 


B = <f> 


SO “ ’'bs 


i 

n 

1 

p 


1/2 1/2 1/2 1/2 
\ - C - f B + C - f 

T /2 T 72 ^ “ 


o a ^^2 1/2 

B - C - I A + C 


a 

2 


c = - V - 1 + ^ 

GS ’'bs p + 4 


a = 


(2esqNA)'^ 


ox 


Finally, the pinch-off contribution is given by: 


4>-(L) 
r n 


(3.4.29) 


Q^(C)d^^ 








P 4 >P 

L&e L 


P(V33-#sl) -P^(L) 

e -e 


2 


t‘ 

e’^ dt 




Thus (3.4.23) , along with (3.4.26, 28, 29) given the drain current 
when the source is in A and the drain pinched off (C). 

If the drain is not pinched off, then the B integral goes all the way 


to the drain, i.e. to %(L) or ^ integral is omitted. 

Similarly if the source is in region B, then the B integral starts at 
♦^(0) or * 20 , and the A integral is omitted. 

Evidently this means that a single current expression is not obtainable 
by this approach. However, in each portion of the device, the equation giving 
the contribution of current is simple. 

The contribution of current from region B of the device, is precisely 
that of Van der Wiele. It should be noted that the approximation we used to 
obtain this result breaks down when the device is strongly inverted (region 
A). It thus seems that a 'best' simple approximation in region A will not 
yield the Van der Wiele result. It is our belief that the treatment of region 
A in this report should turn out to be more accurate, but this MUST be tested 
numerically. A numerical testing of the various expressions for will have 
to be a very careful one, since in strong inversion the boundary layer at x =0 
caused by the high electron concentration there will create problems for a 
general algorithm. 

In conclusion, the strengths of this approach are that different 
approximations for are used through the device, depending on the level of 
inversion. The inclusion of a pinch-off term is also a much needed addition. 

To fully complete this model, one of the following is required. Either 
a numerical testing which shows the value of to be accurate in region A, 
or a method of finding an integrable approximation for x^. , the inversion 
layer thickness. 
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IV, SOURCE AND DRAIN MODELLING 
4.1. Purpose and Prel iminaries 

In the explicit formulae for drain current used by previous clinics for 
parameter extraction, the contribution of source and drain regions was limited 
to the inclusion of (the potential at the source, usually grounded) and 

(the potential at the drain) in the derivation of and the 

potentials at the source-end and drain-end of the channel respectively. 
Theoretical and experimental evidence suggest that this approximation is 
acceptable for long channel devices, but that it is one of the sources of 
disagreement between predicted and measured current in small devices, 
particularly when they are operated at low gate voltage (Vg^ s 2) and 
negative substrate bias (Vg^ S -2). 

This year's clinic has contributed a more accurate analysis of source and 
drain regions and it has at least in part achieved the goal of incorporating 
in the drain current expression two quantities which were previously 
neglected: 

- the source and drain doping profile; 

- the ratio between the depth of the channel and the depth of the source 
and drain regions. 

This result was achieved by applying variational techniques to determine 
the potential close to the ohmic contacts and methods of complex analysis to 
obtain the so called crowding resistance. We shall present here the source 
region analysis. The computations for the drain region are similar and are 
therefore omitted. 


Metal Electrode 
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Q 


Figure 4.1 .1 

Figure 4.1.1 displays the two regions into which the source is divided. 
At the boundary PQ between the two regions the flow lines are assumed to be 
perpendicular to PQ. Region I is the ohmic contact area of the source and 
Region II is the crowding area. We shall present an analysis of Region I 
first. 



Figure 4.2.1 


Figure 4.2.1 shows the partial differential equation and the boundary 
conditions satisfied by the potential $ in region I. We see that at the 
ohmic contact [-r;0] we have a Robin type boundary, while at the insulating 
boundaries the conditions are of Neumann type. The problem suggests that the 
potential t minimizes the functional 
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(4.2.1) T(4>) = 



^ e 


llV4'||^dxdy + X 

1 

«• oo 

J 

0 


+ 00 


4>^(0,y)dy 


with ^=0 outside of [-r^ 0]. 

A first attempt at minimizing (4.3.1) can be tried by assuming that * 
is x-independent . We then obtain 


+00 


(4.2.2) 


(e<i>^ + X«l>^)dy = T(4>) 


Therefore, using Euler's equation, we obtain 

e*yy - X$ = 0 

or 

(4.2.3) vyy 

' ^ = Me + Ne in (-r,0) 

‘ 4> = ay + b in (0» + *) 

^ ^ = cy + d in (-“» 0) 

The constants M,N,a,b,.c.,d, need be selected so that the boundary conditions 
and the continuity requireiiienrs are satisfied. We obtain 

A cosh (7 ^ y+B) 

■ k6y + b 
d 


(4.2.4) #.(y) = 


t 


in (-r,0) 
in (0, + *) 
in (-",0) 


The four constants A, B, b, d satisfy the linear system. 


(4.2.5) 


f A cosh B = b 


— A sinh B = k6 


A cosh (- /- r + B) = d 


e A sinh ( 


/X 

-V t r + 


B) = 0 


Therefore 

(4.2.6) 



This approximation is obviously not optimal. 

A better approximation to 4> can be obtained by assuming 


(4.2.7) ♦(x,y) = m(y) + 4 n(y)(x-e)^ 

t 


Then m and n satisfy the systems of differential equations 


(4.2.8) 


m" + 4 n" = — (m+n) 

— m" +4 n" = n + — (m+n) 

o 0 ^ k £ 
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for X f 0 (in (-r,0)) , and 


(4.2.9) 


m" + ^ n" = 0 


m“ + I n" = -i. n 

5 3e^ 


for X = 0 (outside (-r,0)). Setting 

“-/Fy s = J^r 


we obtain 


%u ' 3 "uu ' " 


(4.2.10) < j 


1 


5 %u = TTx " * " 

in (-S, 0) and 


(4.2.11) 


m +4 n = 0 
uu 3 uu 


%u g ^uu " eX 


outside (-s,0) . Taking into account the boundary conditions on we have 


m = 


4 B,a \ ®3 


3 “2 


n = B.,e 


/15 

VTx " 


(4.2.12) 
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in (0, + «), and 


(4.2.13) 





+ C 


3 




in - s). 

In (-S, 0) the system (4.2.10) gives the following algebraic equation 
for n 


(4.2.14) 6^ - (6+^■) 6^ + = 0 

e A, 6 A. 

For \e < < 1 we have 


6 




+ . . . 


6 


-1 


63 - 1 
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Accordingly (e = ) 

(4.2.15) 

n(u) A-|e®^ + A 26 + A 3 e^ + A 4 e"'^ 




l^m(u) =s -i A^e®^ - 1 A£e - ^^(A 3 e'^ + A^e"^) 


0 u . 2 


I I 

Continuity of m , m , n , n at u = -s and u = 0 gives the following 
system of 8 linear equations in A^, A^, A 3 , A^, B 2 .B 3 , C 2 , C 3 . 

Ai + A 2 + A 3 + A 4 = 83 


IAi + i Ap + ^ A3 + ^ A4 = i Bo - Bo 

0(Ai-A2) “ A^ = "^^2 

ie(Ai-A2) (A3-A4) • -ieBj-kj/f 

-es es -s s 
A^e + A20 + A3e + A4e = C 2 

-es -es _s s 
eCA^e - ) + k^Q - A4e = 0C2 

ie(Aie - Aje ) (Aje - A4e ) = ‘ec 


With a = i 1 
3 e\ 


A 4 = 


2s 

a(e - 1 ) 


A 3 = A 4 e 


^2 = -A 46 


2 s 

S(l- 0) 


3^2 • 

b = k 5 we find that 

B 3 = aA 4 (l+e ) 

Bp - A2+A2'*‘A3+A4 

- 0 s es 
Cp - ” Ape 


(u= 0 ) 


(u=-s) 


2 s s 

A, = k. (l-e)-e ( 1 -t-e) ^3 “ ^ 4 ® 

20 
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Hence 

^1 ~ ^ T kSeVeX 


^2 ~ \ ® 




2 kSsV^sX 


82 


kSsVeX 


k6’ 



^2 ~ ^ k6e\/eX 




s 


+ e 


-s 



Y k6e\/eX e 


-2s 



Consequently 


if y<-r 


.2.16) •l>(x,y) ~ j 


A ^^4 


^ y 

V e 


(12 f - 6(x-e)^) t e 


Vl5 f 

£ 


(3(x-e)^-e^) 


if- r<y <0 


kS 

12 


’-fx 


-V^l T 


T ■*■ + ev4x e -3 


(x-e)‘ 


if 0<y 


The neglected terms are exponentially small. 
The current outflow is (y = 0) 
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(t 


(4.2.17) 


4>ydx = k5e 


^0 


The current inflow is (x = 0, ye[-r,0] ) 
0 


(4.2.18) 


= k6e “ +... 


-r 


and we see that the two flows agree in first approximation. 


The current inflow is concentrated in a boundary layer near y = 0. There 
are two layers, whose heights are and whose widths are and ^ respective 
(see figure 4.2.3) . The solution is independent of r and the current 
enters (or leaves at the drain end) at the corner where the ohmic contact is 
closer to the channel. Hence we can assume that r = « . The current lines 
are represented in figure (4.2.4). 

We can now compute the resistance between the equipotential line ($=0) 

AB (see figure 4.2.1) and the equipotential line PQ. Without going into the 
technical details (see [8]) of the calculations we shall simply mention that 
the resistance is given by 


(4.2.19) 




dc 


ly 


where a is the silicon conductivity, p. the contact resistivity, and 
m c 

where the inner integral is taken along each current line, t|> = c, from 
the region PQ to the contact, (t is the harmonic conjugate of the 
potential $.) 

We then have $(PQ) - $(AB) = R^I and since $(AB) = 0 we obtain 



Figure 4.2.3 



Figure 4.2.4 
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(4.2.20) $(PQ) = Rjl 
4.3 The Crowding Resistance 

As stated previously, this section will be concerned with the analysis of 
region II (see Figure 4.2.2) for the purpose of finding the crowding 
resistance, namely the resistance due to the distortion of the current lines 
near the junction where the source meets the channel of the device [9]. 

Here we will make use of the Schwarz-Christoffel transformation , [12]; a 
conformal mapping technique which maps the interior of a polygon in one plane 
into the upper half of a second plane bounded by the real axis and an infinite 
arc, or vice versa. 



Figure 4.3.1 shows the geometry of the region we are interested in 
transforming with the Schwarz-Christoffel theorem; Figure 4.3.2, instead, 
contains the plane of the straight line into which the boundary of Figure 
4.3.1 has been opened out: 
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2 ,- plane 


Figure 4.3.2 


The differential equation; 
(4.3.1) ^ -c - ^ 


cz 


1 




ca 




is the expression given by Schwarz and Christoffel for the ratio between the 
two infinitesimal vectors dz and dz^. 

Letting z^ * Re^“ , the two constants c and a can be obtained through the 
evaluation of the following integral: 


-h 


(4.3.2) 


f 


* 




/„2 21(1 2Y= 


dz = ic 

1 



s ±ic 

•^+h 

J 

0 

R--0 


f Tt 


ada 


and 


(4.3.3) 


-k 

— 

TT 







fo2 2ia 2 I 

% 


dz = ic 



Re -a 

da 

= ±ic 



lR^e^^“-l - 




^ +k 

_ J 

0 


I 

R-V+eo 


Tf 


da 


0 
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Substituting into (4.3.1) the expressions a= ^ and c = 4^ obtained from 

K ITT 

(4.3.2) and (4.3.3) and integrating, we get; 


(4.3.4) 


k tan 


I 



+ h tan 



Set W - U+iV (with U the potential and V the current flux vector) and W = 
InZj or = e ; then clearly Figure 4.3.1 the current lines follow the 
right and left boundaries of the strip from top to bottom giving a total 
current I = , s being the resistivity of the strip. 

W 

Substituting = e and observing that the line of flow along the y - axis 
is V = I , we find: 




r- — T 

(4.3.5) y = 1 - 

k tanh”^ 


1 


— 


h tanh 


-1 



By equating (4.3.5) to zero, we obtain that the potential at (x=0, y=0) 

must be approximately equal to Xa when a << 1 . The value X= 0.308705435 
has been computed for the parameter X using Newton's method to solve the 
equation: 

(4.3.6) tanh (X^+l)^=-4— , 

(xSi)'^* 


derived by equating (4.3.5) to zero with = ln\a. 

The exact value for U is obtained as an inverse function by setting 


X = 


1T2U 

*.21 



2 2l 



. Then 

=lln 
2 : 

X +a 

J 

J-x^_ 


The substitution of x into (4.3.5) with y=0 plus a little algebra gives 
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the equation: 


1-x ^ ( x-a\ 


1+x 

that we solved using Newton's method with values for a in the 
(0.01,0.034). 

Table 4.3.1 shows the numerical results computed for 
(4.3.7) Uq = InXa and 


(4.3.8) = iln 


Li - 


given the values of X, a and x. 

Now let Uj = U(+<») . With this assumption, equation (4.3.5 
transformed into: 


(4.3.9) y, .2 


2'-'i 2. / ^^1 

k tank"Y V h tanh'^ 1 a ^ 


2U 


1 


2e W 1 


2Uj 2 

2e ^ + a^ 


along the wide part of the strip, or 

2U, 


(4.3.10) 


k fn ^^“ 5 — h fn 
l-a'^ 


2kU 


i.i 


TT n 


k en - h in 


k2-l,2 


k-h 


2kU, 


+ A 


where A = — 


k In -y - y - h «n ^ 


Hence 


^1 ^ 2k * 


interval 


can be 


(4.3.11) 




(approx, values) (approx, values) (exact values) (exact values) 
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Finally, the total resistance of the wide part of the strip in Figure 4.3.1 
not connected to the narrow one is given by: 


(4.3.12) 
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) 


and, since region II of Figure 4.2.1 is exactly one half of the region in 

Figure 4.3.1, (4.3.12) gives the total resistance also for region II, with 

the resistance of only its wide part and -A - u„ the crowding 
^ 2k ° 

resistance. 

We find: 

(4.3.13) * Crowding resistance * 


In particular for a<<l we get: 


1 

2kTT 


k 


4h^ 

k^-h^ 


h lr\ 


k+h 

k-h 


(4.3.14) 


I 

n 


4n2 + ^n 


k 


+ in 




Table (4.3.1) is a comparison between exact and approximate values of both 
and Uq. Figure (4.3.2) illustrates the relationship between Rq , Uq and a 
for the exact and approximate cases. 

4. Remarks 

The procedure described in 4.2 , 4.3 allows us to compute the 
resistance, R^ , of the source region and to express the surface potential, 
♦50 , at the source end of the channel as product of R^ and I . In doing 


so we obviously make an approximation, which is valid at least as long as the 
charge-sheet model holds. We therefore obtain 

(4.4.1) *30 * “s' • 

A similar procedure at the drain gives 

(4.42) Yds - ♦SL = RqI 

where Rq is the resistance of the drain region. 

We can now eliminate and $5^ the explicit formula of Iq and 

we obtain an equation which defines Iq implicitly. 

(4.4.3) Ip = G(IoRs.IpRp.Vss,VBs ) . 

This may seem a big step forward at first glance, but we would like to recall 
that R5, Rq contain the crowding resistance based on the device geometry and 
specifically on k , the ratio between the depth of the channel and the depth 
of the source and drain regions. Therefore k will usually be different at 
the two ends and k=0 at the drain end when the pinch-off regime occurs. 
Thus the surface potential plays an important role in determining the value of 
k. Nevertheless it remains the fact that the doping concentration of the 
source (drain) region and the crowding at the two ends of the channel have 
made their first appearance in the drain current expression. 
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V. NUMERICAL RESULTS 


In sections 5.1 to 5.4, we will describe the numerical algorithms which 
the clinic used for the evaluation of surface potentials (<|)^q and 
in the various models, followed by a numerical comparison, noting accuracy 
and efficiency. 

Section 5.5 gives some notes on the parameter extraction procedures. 
5.1 Exact Surface Potential 

It should be pointed out that all the models mentioned in this report 
calculate and from the same basic one-dimensional model. The 

'exact solution' is the one which uses the full 1-D equation in evaluating 
the surface potentials. The others use approximations to it. 

The equations to be solved are those of ch. 2.2, (see 2.2.10) namely 


5.1.1 <|>^ - Vgg + aF(4.g) = 0 


with F given by 
5.1.2 F(x) = 1 


P(x-Vpc“<l>„) P(4 *p’*’Vdc~^) 

e BS n F BS 


p(|)p P4ip 


e -e 


-e 


and 

5.1.3 



2qn,E3 


1 

2 


To obtain (|)2 q we solve 5.1.1 with (j>^ = (Ji^q and ~ 'l>p " '^bS' 
For we use (|)^ = and (|)^ = + <(>p - Vg^. 

Newton's Method was used for the solution. 
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The Newton parts of the algorithms for and are the same, 

so only the general case will be explained with = x. 

The function whose zero we seek is 


5.1.4 


6(x) = X - + a F(x) . 


The derivative is given by 


5.1.5 


G’ (x) = 1 + 

^ 2F(x) 


with H given by: 
5.1.6 


p(x-V„--<() ) P((|)p+Vdo-x) P4»p -p(|)p 

H(x) = pa(e " -e ^ +e *^-3 


Newton's algorithm, then, becomes: 

G(x^) _ 2F(x^)G(x^J 

5.1.7 ""k+l " ■ G'TV " ’ 2 F(x^)+H(x^) • 


5.1.8 


The Fortran Statement Functions actually used are as follows: 

P(X“Vq<.“iJ) ) P('j*p‘*'^'Dc“X ) P<l»p “p<|>p P<j*r 

FSQ(x) = e ^ +e +p(x-Vg3)(e *^-0 ^-e "^-e " 


5.1.9 


F(x) = DSQRT(FSQ(x)) 


5.1.10 

5.1.11 


G(x) = X - Vg + a F(x) 

P(x-Vb--(|) +tna) p((j)p+\/B--x+ilna) p((j)p+2,na) p(-(|>p+«,na) 

H(x) = p(e " -B ^ +e ^ -Q ^ ) 


In the function H, the tna term in each exponential is equivalent 
to multiplying the expression by a. By using this form, however, the size 
of the exponents is reduced, minimizing the risk of overflow. 
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The general form of the curve G(x) is given below in Figure 5.1.1. 



The key to a successful algorithm based on Newton's Method is the initial 
guess for the zero. Which criteria should we use? 

It was found that Newton's Method could be made to always converge 
for the parameter ranges given in (6.3) under the following scheme. Suppose 
we are given some starting value, x. Then: 

5.1.12 i) Increment x by 0.25 until FSQ(x) > 0 
ii) Increment x by 0.25 until 6(x) > 0 
iii) Use Newton's Method. 


The starting value for was given as 0. (|) 2 g is generally positive 

but can be less than zero. Thus an initial value of 0 will normally converge 
faster since we do not need to step up through all the negative values. 

The starting value for was given as <|) 2 q. Although the approxi- 

mation (^^ 1 ^ = (jjgQ + Vpg might seem natural, it is not practical. When 
the drain is saturated, the approximation is very bad, and the starting 
value will be much too high. In this case, it is observed that Newton's 
algorithm will often converge to the wrong root - the left hand zero of g. 
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If the bounds on the parameter are changed substantially, then the 
surface potential calculations should be checked over the new range. The 
easiest way to do this is probably to make sure that G'(x) is positive 
at the root. 

If^the algorithm does fail under different parameter ranges, there 
are a few things which can be changed before the need to panic becomes 
apparent! 

Firstly, convergence of (|)^q is better guaranteed if the initial guess 
is set to Vgg instead of 0. (With the parameter bounds given in this 
report (5.5) convergence is guaranteed.) 

With this change should also come the following: instead of incrementing 

X until 6(x) > 0, do it until G(x) > 0 and G'(x) > 0. Then it will 
not be possible to start close to the wrong root. 

Another possible reason for convergence to the wrong (left hand) root 
is if the increment in x is too large. Thus use a smaller increment in 
steps i) and ii) of the general algorithm. 

The checks and changes mentioned above are not included in the para- 
meter extraction routines as they stand, since the algorithm always converges 
for our physical bounds on the parameters. Each sophistication costs time, 
and since the programs take a good deal of that commodity already, we have 
kept 'extras' to a minimum. 

5.2 Approximate Surface Potential (Brews) 

The Brews Approximations amount to the following: 
replace the functions in 5.1.8 and 5.1.11 by 
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5.2.1 FSQ(x) 


p(X-Vp--(j5 ) P({)p 

X) = e n ^ P(x-Vbs)® 


P(x-Vp<.-(^ ) 3<j>p 

5.2.2 H(x) = p(e + pe 


Then define the functions 5.1.9, 5.1.10, and the algorithm 5.1.12 is 
seen to converge. 

There is little point in using this algorithm for Brews' surface poten 
tials in a program, however, since the precise solution requires virtually 
no extra effort and convergence is the same. 

A comparison of Brews' results to the precise solution is given in 

5.4. 


5.3 Approximate Surface Potential (Wiele) 

Van de Wiele's solutions for surface potential at source and drain 
ends of the channel are extremely simple to calculate but unfortunately, 
not very accurate. A crucial first step is the determination of the in- 
version level associated to a particular surface potential. The condition 
for strong inversion is given by the relations 

Vgg » Vj and jointly + Vg^ " 'l>p '^1 


Vj = Vg3 - 2<j.p + ix^ - xCix^+Vg)* - I tn((ix^+Vg)"-ix) 


V 


0 




X = (2s^^)"^^ox 


5.3.1 
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and the threshold voltage, Vj, is given recursively by, 


5.3.2 V.|. ^ = Vpg + 2(})p - |A + ^(4A^+vj.-Vpg-Vg2- J) + | «.n 3 + 


I in(-iA+(iA2+v|-Vpg-Vg5- i)") 


(see [ll]eqns. (14), (13), (17)). 

A simple loop is used to find the threshold voltage for each combination 
of Vq 3 and Vg 3 . 

When the relative error between consecutive iterated values of Mj 
falls below 10 it is assumed to have converged to its limit. If it 
is determined that a condition of strong inversion exists then the Van de 
Wiele strong inversion approximation is used. This approximation is a re- 
cursive one and given by 

5.3.3, ^ *0 * "bS * i * I ■ 

(See [11] eqn (9) ) , 

*6 

Here again, a relative error of less than 10 % between iterations is taken 

as proof of convergence. 

If the conditions for strong inversion are not met, the weak inversion 
approximation is used. This is given simply and nonrecursively by 

5.3.4 = Vg3 + iA^ - . 


(See [11], eqn. (12)). 
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5.4 Comments (numerical comparison) 

As Van de Wiele uses the same formulae for the surface potential at 

the source and drain ends of the channel we may use the same approach for 

calculating either one. The only difference between source and drain end 

of the channel lies in our value for <i> . At the source end we take 4> = 

n n 

whereas at the drain end we take = <i>p ■ + V^g. The Van 

de Wiele model has the advantage of being computationally direct, involving 
limited logical checks and only simple loops. Unfortunately, this simplicity 
is at the expense of accuracy. Two runs were made, each for a different set 

Q 

of parameters. For run 1: \ = 2.312, <<,_ = 4, C = 6.9x10 . For run 2: 

• OX 

-8 

\ = .4042, ct>^ = .35, = 23x10 . ctgQ and were calculated through 

the following range of voltages: 2 < Vgg < 5, Vgg « 0 and -2 > Vgg > -5, 

^ - '^DS " 

Each run compared the potentials calculated by the Brews and the Van de 
Wiele models with the potentials calculated by the exact, Pao-Sah model, for 
each possible set of voltages, and expressed the difference in terms of a 
relative error. For the accurate measurement of error, the height of <t>^ 
above Vgg must be taken into account. Thus the expression for relative 
error becomes: 100 x (EXACT) - (APPROX) )/(* (EXACT) - V.-) • 

In the first run. Brews' approximation yielded an RMS error of .3% and 
.5% for and respectively. At no time did the error rise above .7%/1.6%. 

The Van de Wiele approximation v/as not quite as good yielding an RMS error of 
1.4% and .2% with maximum errors of 5.3%/1.2%. For the second run. Brews showed 
0% RMS error in <j>gg and There was only one instance of non zero error 

which was an insignificant .1%. Van de Wiele returned. 1.1% and 3.9% RMS errors 
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with individual cases as high as 2.4% and 1 . 1 %. Clearly then for the pur- 
poses of parameter extraction, the Brews approximation shows itself to be an 
excellent one with almost insignificant error. The Van de Wiele approximation, 
while not as good, is still reasonable for \ = 2.312. However, for small \ 

(x = .4042) errors become noticeable and its use becomes questionable. 

5.5 Parameter Extractions 

The parameters used by last year's clinic with the Brews model have been 
modified by the clinic. Previously it was the case that 2 of the parameters 
each contained two required quantities. Namely 

5.5.1 P 2 = P2('*’f’^ox^ 

P 5 " ^5^'^'^'-’^ox^ 

Since <|)p, and W/L are the real quantities we wish to find, we have 
set the new parameters to be: 

Pi = 

5.5.2 P 2 = 

P 5 = W/L 

What bounds should be put on the paramters? From physical considerations, 
we have the following; 

5.5.3 0.29 s <jip s 0.41V corresponding to 10^^ 5 5 10^^ > 


-0i95< Vpg s -0.8V 

for 

n-channel , 

-0.3 < Vpg < -0.2V 

for 

p-channel 
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W and L are given for each device. Acceptable bounds for W/L 
would be 10% to either side of the given value. 

For the Brews and Van de Wiele models, the full one-dimensional solution 
for surface potentials is used. Also, we have coded Van de Wiele's approxi- 
mate solutions in the hope that they may be more computationally efficient. 

A comparison of the various surface potential calculations was given in 5.4. 
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6. CONCLUSIONS 

The new formulation of drain current given in (3.4) is, we feel, more 
acceptable than those previously derived. If an approximation of x^. could 
be found in region A, then it would be complete. As it stands, however, the 
importance of such a correction is not known. Tests should be carried out 
comparing the approximations for with numerically computed values — this 
should also be done for the Van de Wiele model. 

What is it about this model which is more acceptable? 

Firstly it includes the contribution of current from the pinched off part 
of the device which is not seen in any of the other models we have looked at. 

1 . 

Secondly, in the Brews and Van de Wiele models, the result Or. = 

is used for the depletion charge. This result is the same as the clinic's 

model in region B, but look at how they calculated this value. Firstly, they 

assume that all the charge not due to electrons is due directly to the doping, 

N^, and that the depleted region is uniform to a depth W, beyond which charge 

2 2 

neutrality exists. Thus they solve the equation d <^>/dx = q N^. In fact, 

as we go down into the bulk of the device, the hole charge, p, increases con- 
tinuously and only asymptotically to the value N^, which gives charge neutrality. 

The boundary conditions they use at the surface are those of the full 1-D Poisson. 
Hence this approach v/ould seem to neglect the sharp increase in the actual 
potential near the surface caused by the inversion layer electron charge. The 
clinic has not been able to justify the use of this approach for finding Qg. 

In our model we allow all the charge distributions to be continuous, and 
do not lose the inversion layer structure. 



Brews claims that his model is valid in all regimes of operation. Numeri- 
cally, this is notwhat we find. By using parameter extraction, we can get a 
good fit to wide ranges of data. However, there are certain combinations of 
voltages which without fail cause large errors. Notably '^GS " '^BS " -3. -4, -5. 

The departure of the model from reality in these conditions can be by 30 
or 40% compared to less than 5% for most all the other voltages (for one given 
fit to data). It would thus seem reasonable to assume that the model is having 
trouble fittingto these particular applied voltages. The reason for this may 
be that his model uses a single formula for whatever the operating condi- 
tions at any point in the device. The clinic, however, splits the device into 
3 regions with different approximations in each. 

What is required now is a thorough testing of the various forms of 0^ 
which are used in the two models, as well as an attempt to find x^. in region 
A as an integrable function. 

All the models mentioned in this report are based on the Pao-Sah formula- 
tion of drain current. Typically, instead of integrating directly to find 
Q^,the route using " *^n taken. 

As has been noted previously, the Pao-Sah double integral formula is found 
to be highly accurate for long-channel MOSFETs. For short-channel devices, 
accuracy is not guaranteed. 

The main inaccuracy in this model is due to its main feature: namely the 

integration of Poisson's equation. In order to do this integration, two assump- 

2 2 2 2 

tions are needed: (i>^ = ■i)p(y) and a 4>/dy « a 4>/bx . The effect of these 
assumptions is to make the model a combination of two uncoupled one-dimensional 
models, rather than a true two-dimensional model. 
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What is needed now is to start with a model based on Pao-Sah, and to 
extend it to something which is "more two-dimensional". Part of this 
approach should be in the source-drain modeling. Thus and are 

no longer calculated from the one-dimensional Poisson, but from full solutions 
to the source and drain equations. The clinic has analysed much of the large 
amount of work done in this area. We have chosen the approaches which are 
most mathematically sound (yet not too complex to be unworkable) and presented 
them as a full and concise reference for future attempts to include source 
and drain effects in a drain current model. 

Along with this approach, we need to add some y dependence to the solu- 
tion of Poisson's equation. This is particularly important at the source and 
drain where irregular device geometry and rapid changes in carrier concentra- 
tions along the channel (e.g., caused by pinch-off) can make the one-dimensional 
approximation very unrepresentative of the true situation. 

A clear understanding of the tricks and techniques used in one-dimensional 
modeling is paramount before attempting the step to quasi -two-dimensionality. 

The main achievement of this Clinic is that we have clarified the approaches 
commonly used in one-dimensional modeling, consolidated much of the large amount 
of work done on source and drain modeling, and possibly produced an optimal one- 
dimensional model for MOSFET current. Initially, we had planned to produce a 
two-dimensional model, but the lack of what we considered to be a fully justi- 
fied 1-D formulation led us, instead, to concentrate on this, thus producing a 
good basis for future mathematical modeling in this area. 

For the computing, we have produced efficient, accurate algorithms for the 
solution of surface potentials. These have been incorporated into parameter 
extraction programs for the Brews and Van de. Wiele, madeTs:. 
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APPENDIX I 
NOTATION AND SYMBOLS 

The notation and symbols detailed below have been used exclusively 
throughout the report. 

L Channel Length 

Lg Effective channel length 

W Width of the device 

x^(y) Depth of the inversion region 

D(y) Depth of the depletion region 

IJ.^ Electron mobility 

Effective electron mobility 

Oxide layer capacitance per unit area 

n^. Intrinsic Carrier Concentration 

N^ Dopant ion concentration (acceptors) 

Charge density due to free electrons in the inversion layer 

Qg Total space charge density 

Q|^ Depletion region charge 

J y component of current density in the channel 

y component of electron current density in the channel 

n Carrier density of electrons 

p Carrier density of holes 

<t>(x,y) The potential measured relative to the source 

<t>g(y) <ji at X = 0, the surface potential 

(|> 2 Q *^*s^*^^ i-e., surface potential at the source 
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potential at the drain 

4>j^ Electron quasi Fermi potential 

<^p Hole quasi Fermi potential 

-1 

= |3 bulk Fermi potential 

6q Permittivity of free space 

<5 Dielectric constant for silicon 



'BS 


'DS 


'GS 

^GB 


'FB 


ox 


Dielectric constant for oxide 

Boltzman constant 
Absolute temperature /®K 
= q/kT 

, . . . . i 

= Ug/ in^.qB))": Intrinsic Debye Length 
= (cg/NAqp))^: Debye Length (bulk) 

' '®'”a'-b'^ox 

body bias relative to the source 
drain bias relative to the source 

gate bias relative to the source less the flat band voltage, Vpg 

gate bias relative to the body less V^g 

Flatband voltage 

Thickness of the oxide layer 

= (°n,Pc) ^ where is contact resistivity, is material 

conductivity 
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APPENDIX II 

NOTATION AND SYMBOLS OF VARIOUS MODELS 


This section is a translation for changing expressions found in some of 
the papers we have used into the clinic's standard notation. It should be 
noted that the source to substrate bias is zero in Pao-Sah [6] and Brews [2]. 
Thus it is only valid to set Vg^ = 0 when comparing these to the other three 
papers. 


Clinic 

Pao-Sah [6] 

Van de Wiele [11] 

Brews [3] 

Brews [2] 

'^BS 

0 

\-h 

~^BS 

0 

'^GS 

u^/p 

h-S- ''fb 

Vq - Vgg or Vgg 

h 

'^DS 

Uq/P 


Vds 

h 

0 

0 

Vc 

0 

0 

4> 

u/p 

0 - Vg 

^ ■ '^BS 

4> 


Ug/P 

^s-Vs 

■ '^BS 

h 


Up/p 

*F 

<i>F 

4>q 

'^n ^ '^BS 

(5+Up)/3 

v(y) + 

“^F^ ■ ^BS 

^f 

X 

X 

X 

y 

X 

y 

y 

y 

X 

y 


— 

*^sc 


<5s 

Od 

— 

% 


— 

'5n 

— 

% 

Nl 
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